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ABSTRACT 

Star formation rates (SFR) larger than 1000 M yr _1 are observed in extreme star 
bursts. This leads to the formation of star clusters with masses > 10 6 Mq in which 
crowding of the pre-stellar cores may lead to a change of the stellar initial mass 
function (IMF). Indeed, the large mass-to- light ratios of ultra-compact dwarf galaxies 
and recent results on globular clusters suggest the IMF to become top-heavy with 
increasing star-forming density. We explore the implications of top-heavy IMFs in these 
very massive and compact systems for the integrated galactic initial mass function 
(IGIMF), which is the galaxy- wide IMF, in dependence of the star- format ion rate of 
galaxies. The resulting IGIMFs can have slopes, 0:3, for stars more massive than about 
1 Mq between 1.5 and the Salpeter slope of 2.3 for an embedded cluster mass function 
(ECMF) slope (/3) of 2.0, but only if the ECMF has no low-mass clusters in galaxies 
with major starbursts. Alternatively, (3 would have to decrease with increasing SFR > 
10 M Q yr _1 such that galaxies with major starbursts have a top-heavy ECMF. The 
resulting IGIMFs are within the range of observationally deduced IMF variations with 
redshift. 

Key words: stars: formation - stars: luminosity function, mass function - galaxies: 
evolution - galaxies: starburst - galaxies: star clusters - galaxies: stellar content 



1 INTRODUCTION: CLUSTERED 
STAR-FORMATION AND THE IMF 

The stellar initial mass function (IMF) is one of the fun- 
damental astrophysical distribution functions. It defines the 
ratio of low-mass stars, which do not contribute to the chem- 
ical evolution over a Hubble time but lock-up baryonic mat- 
ter, to high-mass stars, which power the interstellar medium 
and enrich it with metals through AGB-winds and super- 
novae. It further determines the mass-to-light ratios of stel- 
lar populations and influences the dynamical evolution of 
star clusters and whole galaxies. 

The low-mass end of the IMF is found to be indepen- 
dent of environmental influences like metallicity and den- 
sity, which may be understandable theoretically by a nearly 
constant Jeans-mass due to the way how molecular cooling 
rates scale with density and temperature (jElmcgrcc n et all 
2008). In the regime around and below the hydrogen burn- 
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ing mass-limit it now seems that brown dwarfs form an in- 
divid ual distribution disjoint but related to the low-mass 
IMF (|Thies fc Kroupall2007l , l2008l ). The formation and the- 
oretical basis of the IMF of massive stars (> 10 Mq) is 
less well understood with at least two competing theo- 
ries (competitive accretion vs. single star accretion) having 
been developed (jBonnell et aflll998l. 120041: iBonnell fc Bate! 
120061 : iTan et al.ll2006l : iKrumholz et al.ll2009l ). Observation- 
ally, though, the slope of the high-mass IMF within star 
clusters seems to be as independent from the environment 
as the low-mass slope. Unresolved multiple systems have 
been shown to have no effect on the m > 10 Mq IMF 
|Mafz Apellanij|2008l : IWeidner et"ai] 120091 ) but for m ^ 1 
Mq unresolved multip le stars effect the o bserved slope of 
the IMF substantially (|Kroupa et al.lll99ll ). 

Over the last years it has become clear th at star forma- 
tion takes place mostly in embedded clusters (|Lada fc Ladal 
2003). each cluster containing a dozen to many million of 
stars (Kroupa 2005). Within these clusters stars appear to 
form following the canonical IMF, £(m) oc m~ a , with a slope 
of 1.3 for low mass stars and the Salpeter/Massey-slope of 
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2.35 for massive stars (for more details on the canonical IMF 
see Appendix [5} . 

Not only stars follow a mass function but also (young, 
embedded) star clusters. The embedded cluster mass func- 
tion (ECMF) has been found to be a power-law, £ cc i oc , 
with a rather constant slope of « 2 for largely different envi- 
ronments from the quiescent solar neigh bourhood to the vig- 
orously star-forming Antennae galaxies (|Lada fc Ladall2003l ; 
iHunter et af1l2003l ; lziiang fc Falllll999l) . 

The upper mass end of the ECMF (the most-massive 
young cluster within a star-forming gala xy) has been 
found to depend on the SFR of the galaxy (| Weidner et al.l 
120041 ; iBastianl | 2008). Furthermore, it appears that that 
the mass of the most-massive star in a star cluster is 
related to the mass of the cluster in a non-trivial way 
l|Weidner et al.ll2010h . Low-mass clusters seem to be unable 
to form very massive stars. The physical reason for this 
empirical relation has not yet been found. But it can 
be presumed that the interplay between stellar feedback 
(ionising radiation and winds) from the massive stars and 
the gravitational potential of the star-forming cloud might 
play a role. In this picture, the feedback of the massive 
stars overcomes the binding energy of the cloud and 
star-formation is terminated. Such a process would directly 
couple the mass of the cloud with the feedback of the stars 
which is directly correlated with their mass (jElmegreenl 
Il983h . The majority of these embedde d star clusters then 
dissolve quickly due to gas-expulsion (IKroupa et all 200 ll; 
Adams fc Mversll200ll; IKroupa fc Boily||2002j; lLada fc Ladal 



20031; Ide la Fuente Marcos fc de la Fuente Marcos! |2004 
Parmentier et a l. 2009). 



A direct consequence from clustered star-formation is 
that the composite stellar population in a galaxy, which 
results from many star-forming events, is the sum of the 
dissolving star clusters. Thus the integrated galactic initial 
mass function (IGIMF) is the sum of al l the IMFs of all the 
star clusters ( Kroupa fc Weidner] 120031 ; IWeidner fc Kroupal 
l2005l . see also lvanbeverenlll982f ). 



£iGiMF(m,t) = / £(m ^ m max (M cc i)) 

J Af ocl min 

•£ecl(M ec i)dM e cl, 



(1) 



where M cc i iin i n is the minimal embedded cluster mass, 
M ec i,max(S'7 ? _R(t)) is the maximum embedded cluster mass 
which is dependent on th e SFR of the galaxy and is given 
by the following equation (|Weidner et al.ll2004 ). 



Meo, max = 8.510 4 x SFR 075 M Q . 



(2) 



And £ecl(M ec i) is the ECMF0, with dN c , 



£eci(Afeci) dM cc i oc M ec f dM cc i being the number of just 



1 Note that, strictly, £ oc i describes the distribution of star- 
forming molecular cloud cores containing only the stellar mass 
formed. The "embedded clusters" in eq. [TJ do not, under any 
circumstances, mean bound or radially well-defined stellar en- 
sembles. Rather, eq. [T] is an integral over all locally correlated 
star forming events, a small fraction of which will hatch from 
the clouds as bound clusters, while the majority disperse within 
about 10 Myr. 



formed embedded clusters with stellar mass in the inter- 
val M cc i, M cc i + dM cc i- Finally, £(m < m max (M GC i)) is the 
IMF within an embedded star cluster with the total stellar 



pend ent on M cc \ l|Weidner fc Kroupal 20061; Weidner et al.l 
l2010h . Note that eq. His originally l|Weidner et al.l I2004T ) 
only fitted to SFRs of up to 4 Mq yr _1 and extra-polated 



for starbursts with SFRs > 10 M, 
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1 in this work. But 
that eq. [5] also 



it has been recently shown dBastia: 
holds for higher SFRs. 

The IGIMF is steeper than the individual canonical 
IMFs in the actual clusters, hereby immediately explaining 
why 03, Held = 2.7 > 03 = 2.35, wher e 03. field is the slo pe of 
the IMF derived bv lScald (Il986l ) and lReid et~ai1 (2002) from 
QB star counts in the Milky Way field, and 03 = 2.35 is the 
ISalpeterl (|l955l ) index. This is due to the fact that low-mass 
star clusters are numerous but can not have massive stars. 

The high mass part of the resulting IG IMF is strongly 
depe ndent on the SFR of a galaxy (Haa s fc Anders! 
|2010[ ). This result has been found to be naturally (with- 
out parameter adjustments) ab le to explain the mass- 
metallicity relation of galaxies ijKoppen et al.l [2007), the 
alpha-element abun dances as a function of galaxy mass 
llRecchi et alj [2 009 ) and the Ha cut-off in star-forming 
galaxies l|Pflamm-Altenburg fc Kroupal I2008T ) . It also pre- 
dicted a discre pancy between SFRs de r ived f rom UV- 
and Ha-fluxes (|Pflamm-Altenburg et all 120071 . |2009|) . a 
result which has recently been confi rmed qualitativ ely 
l|Meurer et alj|2009l ) and quantitatively (|Lee et al.ll2009l ) by 
observations. Furthermore, the Ha-SFR relation calculated 
in the IGIMF theory le ads to higher SFRs of Ha faint star- 
forming dwarf galaxies (Pfi amm-Altenburg et al.ll2007l ) and 
reveals a simple linear relatio n between the total neutral gas 
mass and SFRs of galaxies (Pfl amm-Altenburg fc Kroupal 
120091) . 

The IGIMF theory (eq. [TJ is based solely on observed 
correlations (m max vs M cc i, M oc i,max vs SFR) and distribu- 
tion functions (IMF in the correlated star formation events, 
or "embedded clusters", and the mass function of stellar 
masses formed in these events, the ECMF). It is these ob- 
served correlations and distribution functions which contain 
the relevant star-formation physics. It is this basis which al- 
lowed the I GIMF theory to be so su ccessful since its first 
formulation (|Kroupa fc W eidner 2 0031 ) . The empirical basis 
of the IGIMF theory is important because star-formation 
theory is not advanced enough to make reliable statements 
on the galaxy- wide IMF, and would in any case lead directly 
to the IGIMF theory since star-formation theory will have to 
account for the observed correlations and distribution func- 
tions. 

In the case of very high SFRs (> 30 M©/yr) the 
Afeci.max — SFR- relation results in very massive star clusters 
(M ec i > 1O 6 M0) which may be the progenitors of the Ultra 
Compact Dwarf galaxies (UCD), as th ese most-likely formed 
as co mpact objects of a few pc size (|Dabringhausen et al.l 
2009). And indeed recent HST/Spitzer observations found 
very young massive (1O 9 M0) objects with radii below 100 
pc in low-redshift L yman-Break-Analog g alaxies with high 
star-formation rates l|Overzier et al.ll2009h . 

In this contribution we explore the consequences of very 
high galaxy-wide SFRs on the IGIMF by taking into account 
recently acquired evidence for top-heavy IMFs when UCDs 
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and possibly also globular cluster are formed. This advances 
the IGIMF theory into a physical regime of cosmological 
significance si nce very high SFR s are typically observed at 
high redshifts (jTresse et alj|2007f ). 

The observational situation regarding top-heavy IMFs 
and crowding in massive proto-clusters is discussed in Sec- 
tion [2] while the top-heavy IGIMF model is presented in 
Section [3] The following Section [4] then shows the results 
of the model calculations which are discussed in Section [5] 
The canonical IMF within star clusters used throughout this 
work is presented in Appendix 1X1 



2 OBSERVATIONAL SITUATION 
2.1 Top-heavy IMF 

Despite all the evidence for a universal IMF, recently sev- 
eral indications have begun emerging for a possible depen- 
dence of the shape of the IMF on environment. Several ob- 
servational and theoretical indications suggest the IMF to 
become t op-heavy under ext reme starburst conditions. For 
example, Bau eh et al.l <|2005h find that a natter high-mass 
slope of the IMF at higher redshifts is nee ded to explain 
the ob served numbers of sub-mm galaxies. IStanwav et al.l 
(2005) observe rather steep slopes in the spectral energy dis- 
tribution of high-redshift objects in the Hubble Ultra Deep 
Field. They conclude this can be either explained by a 50% 
higher star formation rate than the current calibrations yield 
or a top-heavy IMF at very hig h-redshifts (e.g. z > 10 - 
15). This is actually predicted bv lClarke fc Brom m (2003). 
A similar evolution of the IMF with redshift such that at 
high-z star-formation is biased towards more massive stars 
is found by studying the amplitude of_the galaxy-stellar- 
mass-star-formation-rate relationship (|Davel l2~008h and by 
comparing the rate of the luminosity evolution of massive 
early-type galaxies in cluster s to the rate of their colour 
evolution (Ivan Dokkum|[2008). As star-formation rates are 
also higher at high-z this can be seen as an indication for a 
top-heavy IMF in starbursts. When integrating the redshift 
evolution of the SFR over cosmo logical time-scales w ithin 
the A cold dark matter cosmology I Wilkins et al.l (|2008l ) find 
that they need a top-heavy IMF with ay = 2.15 to re- 
produce the present-day stellar mass density. However, in 
order to explain the stellar mass density at higher red- 
shifts an even flatter slope is needed. Other constrains also 
come from the element abundances observed in galaxies. 
iFinoguenov et al.l (|2003l ) study the element abundances in 
groups and clusters of galaxies with the use of XMM-Newton 
and ASCA. They show "that while the metal production in 
groups could be described by a stellar population with a stan- 
dard local initial mass function, clusters of galaxies require 
a more top-heavy IMF. " From modelling the iron evolu- 
tion in galaxy clusters and the field a nd co mparing the re- 
sults with observations, iLoewensteir] (2006) concludes that 
a top-heavy IMF is indeed needed in galaxy clust ers but 
not for galaxies in the field. iTornatore et al. (|2004T ) find a 
similar result in their A^body-Tree + SPH simulations of 
a galaxy cluster. A Salpeter IMF produces subsolar [O/Fe] 
abundances. They can explain solar [O/Fe] if the IMF is 
top-heavy for z > 2 and Salpeter afterwards. To account 
for the [q/Fel-cr-rel a tion o bserved in local elliptical galax- 
ies ICalura fc Mend (|2009l ) also require an top-heavy IMF. 



iBallero et ail (|2007al lbh infer the need of a top-heavy IMF 
in the Galactic bulge in order to explain the difference in 
the [O/Mg] to [Mg/H] ratios in bulge Gi a nts co mpared to 
the Galactic disk. But iMc William et all (120081 ) find that 
metallicity-dependent stellar yields for massive stars equally 
explain the ratios. Hints for a top-heavy IMF are also found 
in the centre of the Milky Way. By examin ing the evolution 
of dis k stars orbiting a central black hole lAlexander et al.l 
(2007) found strong evidence for a top-heavy IMF in or- 
der to explain the observed ring of massive stars orbit- 
ing about O.lpc around the Galacti c centre. Using Chan- 
dra X -ray data of the Galactic centre iNavakshin fc Sunvaevl 
(2005) found 10 times less X-ray emission from low-mass 
stars than expected if t he observed mas sive stars formed 
from the canonical IMF. lElmegreenl j|2005h did an extensive 
literature study of Galactic and extragalactic observations 
and concluded that dense star-forming regions like star- 
bursts might have a slightly shallow er IMF, a view share d by 
lEisenhauerl (|2001al Jbh. Furthermore. lBrown et~aH (|2007n find 
that the majority of observed hypervelocity stars ejected by 
the super-massive black hole in the centre of the Galaxy have 
a mass around 3 or 4 M© whilst lower mass stars are lack- 
ing from observations. The high number of carbon-enhanced, 
s-process-enriched unevolved stars among extremely metal- 
poor stars in the halo of the MW has been interpreted by 
iLucatello et~aH l|2005h as evidence for a top-heavy IMF dur- 
ing the formation of the halo. A recent study found that 
the central regions of disturbed galaxies host far more su- 
pernovae lb and Ic (originating from very massive stars) 
while in all other enviro nments supernova typ II dominate 
l|Habergham et alllioioh . This can be attributed to a top- 
heavy IMF in starbursts after galaxy interactions and merg- 
ers. 

Marks, Kroupa & Dabringhausen (in preparation) anal- 
yse the stellar content of globular clusters and the dynamical 
M/L ratios of UCDs using A-body models and deduce that 
the IMF appears to become top-heavy with increasing den- 
sity of the forming object. Given the mass-radius relation 
inferred from the proto-globular clusters and proto-UCDs 
this translates into a dependence of the IMF slope Q3 on 
the object's mass. The application of the IG IMF theory to 
high g alaxy- wide SFRs is based on this result. [Papadopoulosl 
(|2010T ) shows that the molecular gas is likely to have a higher 
temperature in massive star bursts as a result of heating 
through supernova generated cosmic rays. This may be a 
theoretical reason for expecting IMFs that are top-heavy in 
such systems. 

2.2 Crowding in very massive pre-cluster cloud 
cores 

A substantial fraction of stars seems to form i n star cluster- 
like c orrelated events (lAdams fc M yers 200 ]]; lLada fc Ladal 



their final stages of collapse (r ec 


1 S 1 pc. iTesti et al.l 


1998; 


Kroupa 20051; Gutermuth et al. 


2005; Testor et al. 


2005; 


Rathborne et al.l 2009; Schcepmakcr et al. 20091) and are 



therefore in very dense configurations. 

Also, strikingly similar values are found for the size 
of pre-stellar cloud cores, despite large differences in mass. 
iFuruva et all (|2006l ) find radii o f about 5000 AU for a low- 
mass star-forming environment. IfBacmann et all (|200lh find 
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5000 to 10000 AU and Ivan der Takl (|200Ch find 20000 AU 
for massive star progenitors. 

Furthermore, within sta rbursts extreme star forma- 
tion r ates (> 100 MjJyr" 1 , lYan et all 12004 iDaddi et all 
12004 IChqi et alj 12009) are reached, e.q. in the case of 



Ar p 220 (IWilson et ai]|2009\ The recent Herschel study 
by iMagnelli et all |201C| ) of submillimetre galaxies (SMG) 
found that for SMGs with fluxes S(850^im) > 5 mjy the 
median SFR is 960 M yr" 1 . The highest 5(850Atm) flux 
in their sample is about 10 mjy, corresponding to a SFR of 
several 10 3 M yr -1 . 

This leads to the formation of very massive star clus- 
ters according to eq. [2] For the characteristic time-scale 
of conversion o f the interstellar medium to stars, St, of 
about 10 Myr (|Egusa et alj |2004 l2009h , the total mass 
in stars is Mtot = SFR x St. But as also Mtot = 
J tf " | '" m " M «:i(cci(M l; d)(iM ( , c i 1 M cc i imax increases as well 

when the SFR increases. The normalisation of £ e ci(M ec i) is 
obtained by stating there is one most-massive cluster, 1 = 
Jm i Coci(M ec i)dM ec i, forming during each 10 Myr inter- 
val. 

Because young pre-star clusters are compact (r ~ 1 
pc) and since the contracting pre-stellar cores within them 
have dimensions of a few thousand AU the pre-stellar 
cores are likely to be gin interacting and c oalescing above 
a certain mass limit. iBonnell et alj (|l998l ) predict an in- 
fluence on the IMF of the for ming star cluster, and so do 
lElmegreen fc Shadmehrl (|2003h who investigated this ques- 
tion of crowding further. They concluded that crowding 
should be relevant in the most massive clusters (e.g. pro- 
genitors of globular clusters) leading to a top-heavy IMF, 
because the crowded medium would have a higher low-mass 
limit for star formation and therefore more m assive stars 
need t o be build to result in the same total mass. lElmegreenl 
l|2004 ) constr ucted a three-com ponent IMF model from 
these results. IShadmehril (|2004 ) further showed that the 
effect of crowding must have an influence on the IMF in 
starburst clusters. While lElmegreenl |2004 ) and [Shadmchri 
l|2004 ) present models for the IMF in a starburst they do 
not discuss the overall integrated IMF for a whole galaxy 
with a starburst. Certainly, in massive clusters, or UCDs, 
the IMF will be effected by crowding (higher M/L ratio, 
Dabringh ausen et alj|200"9l) . but the question remains if this 
change is large enough to influence the IGIMF significantly. 

In Fig.[T]the crowding is illustrated. It shows how many 
proto-stars of a fixed pre-stellar cloud core size can be fitted 
into a spherical proto-cluster with a radius of 1 pc. This is 
done by taking a range of proto-cluster masses and divid- 
ing these masses by the mean mass of the canonical IMF 
(m m ean « 0.36 Mq) in order to get the expected number of 
stars, -/Vexpec, in the cluster. N eyLpec is then multiplied by the 
volume of three different pre-stellar cloud core sizes - 1000 
AU (thick solid line), 5000 AU (thick short-dashed line), and 
10000 AU (thick long-dashed line). Fig. [T]shows the volume 
of the proto-cluster divided by the total volume of the pre- 
stellar cloud cores. For clusters with more than « 10 5 M 
crowding will be a problem. It is therefore appropriate to 
assume that such high-density environments may have an 
influence on the IMF within such star clusters. 

An alternative or additional mechanism which may in- 
duce star bursts leading to top-heavy IMFs may be through 
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Figure 1. The volume of an embedded star cluster (V oc i, with 
a radius of 1 pc) divided by the volume of all pre-stellar cores 
within it (N s tars * V^tar) versus the embedded cluster mass 
(Af cc i). The thick solid line indicates the assumption of a pre- 
stellar core radius of 1000 AU, while the thick short-dashed line 
assumes a pre-stellar core radius of 5000 AU and in the case of 
the thick long-dashed line, a pre-stellar core radius of 10000 AU 
is used. The value when the volume of all stars equals the volume 
of the cluster is indicated by a horizontal solid line. 



the heating of the molecular ga s by supernova generated 
cosmic rays l|Papadopoulosll2010l ). 



3 THE MODEL 

While the above discussion suggests that the IMF may be- 
come top-heavy, it remains, unfortunately, unclear as to how 
a systematic v ariation of the IMF with increasing M ec i ou ght 
to be realised. lElmegre"enl (|2004l ) and lShadmehrU (|2004 l al- 
ready presented descriptions of the change of an IMF in 
starbursts. They based their models on a log-normal IMF 
l|Scalolll986l ) but we use the multi-power law formulation of 
the stellar IMF. 

The modelling of the IGIMF starts from a SFR which 
determines the upper mass limit of the ECMF according to 
eq. [2] In order to calculate the total amount of stars formed 
in a star-forming period the SFR is multiplied by a time- 
step which is assumed to be 10 Myr. With the use of the 
total mass the ECMF can be normalised for the time-step. 
The ECMF is then divided in 1000 logarithmic bins. As 
the embedded cluster mass limits the most-massive star in 
each cluster, the upper limit of the IMF in each ECMF- 
bin can be calculated. As the number of clusters per bin is 
given by the normalised ECMF, the IMFs of each ECMF- 
bin can be multiplied by the number of clusters in order to 
give the integrated galactic IMF of the time-step. The slope 
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of this IGIMF above 1.3 Mq can then be computed by a 
least-squares fit to the calculated 1G1MF. The value of the 
so-derived a-j, gives a rough indication of the top-heaviness, 
while detailed astrophysical parameters would be derived 
from using the actually computed IGIMF. 

In order to include the possible effect of crowding into 
the framework of the IGIMF theory (§ [2T2J) the results of 
Marks, Kroupa & Dabringhausen (in preparation) on the 
high-mass IMF slope in globular clusters and UCDs are 
used. They derive the following dependence of the IMF slope 
for stars more massive than 1 Mq and for clusters with ini- 
tial masses M ccl > 2xl0 5 Mq, 

a 3 (M ccl ) = -1.67 x log 10 (j^g^ + 1-05. (3) 

The limit of 2x10 Mq is chosen because for clusters with 
masses below this limit Marks, Kroupa & Dabringhausen (in 
preparation) results infer a Salpeter slope. Note that while 
the parametrisation in eq.[3]rests on empirical globular clus- 
ter data which have been corrected for stellar and dynamical 
evolution, the parametrisation we adopt in eq. [3] is not to 
be seen as an established dependence of ag on Af ec i- Rather, 
while giving a good clue as to how the IGIMF may change 
with increasing SFR, the calculations performed here are 
independent of the parametrisation and may thus be easily 
adopted for different parametrisations. 
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Figure 2. Three different IGIMFs from Tab. [TJtogether with the 
input canonical IMF (solid line, eg, IA2j l, All three models have f) 
= 2.0 and for the dotted line SFR = 10 AfQyr - 1 and M ec i m ; n = 
1000 Mq, for the short-dashed line SFR = 100 M yr _1 , M cc i, min 
= 10000 Mq and for the long-dashed line SFR = 1000 M yr- 1 , 
Mecl.min = 100000 Mq. All IMFs are normalised to the same 
total mass in stars (10 11 Mq). 



3.1 Model parameters 

Besides the change in the IMF for very massive clusters, the 
most important parameters in our model are the SFRs of 
the models (which give the upper mass limit of the form- 
ing clusters, M cc i jma x), the slope (P) of the ECMF, and the 
lower mass limit of the ECMF, Meci.min- For this study of 
starbursts four SFRs (10, 100, 1000 and 10000 Mq yr -1 ) 
are chosen. Below 10 Mq yr -1 little impact on the IGIMF 
from the top-heavy IMFs in very-massive clusters is to be 
expected as only very few or none of them are formed ac- 
cording to eq. [2] As already discussed in the introduction, 
SFRs of more than 10 00 Mq yr -1 have recently been found 
l|Magnelli et all l201Ch . The upper limit for the SFRs in 
the models is therefore chosen to be 10000 Mq yr - . The 
value of the slope of the ECMF is usually given with a j3 
around 2 with errors between 0.2 and 0.5 ijLarsenl [2009h . 
thou gh some studies fi nd systematically flatter slopes like 
1.8 jPowell et all 120081 ). And the mass spectrum of giant 
molecular cl ouds, the precurs ors of star clusters, shows a 
slope of 1.7 (|Rosolowskvll2005l ). As the identification of em- 
bedded clusters in extra-galactic objects is very challenging, 
we allow a greater range for this parameter and study f3 = 
0.5 to 2.35. The lower mass limit of the ECMF, M oc i, min , 
is even worse constrained than the slope of the ECMF. We 
therefore study a broad range of possible M oc i, m j n from 5 to 
10 5 Mq. 

The lower mass limit of the canonical IMF in each clus- 
ter is fixed to 0.01 Mq an d the upper mass limit is deter- 
mined by the cluster mass jWeidner fc Kroupalliooi ). 



4 RESULTS 

Tab. [1] shows the IGIMF slopes for stars more massive than 
1.3 Mq for different star-formation rates, different values 
of the ECMF slope, /3, and for several lower limits of the 
ECMF, M ec i,miii- For clusters above 2xl0 5 Mq the IMF 
slope for stars above 1 Mq, 03, is changed according to eq.[3] 
while below 2x10 Mq 03 is kept constant at the canonical 
(Salpeter) index 2.35. For «3 a lower limit of 1 is used be- 
cause clusters with an 03 of 1 or less contain more than 90% 
of the mass in stars more massive than 8 Mq. When these 
stars explode as supernovae the resulting extreme mass-loss 
would completely disperse the star clusters. In such a case 
no globular clusters would survive until today. For clusters 
more massive than ea 1.1 xlO 6 Mq this limit is reached and 
a?3 is fixed to 1.0 for these objects. 

In Fig.[5]three IGIMF examples from Tab.Q]are shown 
together with the input canonical IMF. Depending on the 
parameters a strong top-heaviness can be seen. 

In Fig. [3] the results of Tab. [1] are summarised for ft 
= 2.35, 2.0 and 1.6. Generally, the IGIMF slope decreases 
with increasing SFR for a given /3. This is due to the rela- 
tion between the maximal cluster mass and the SFR (eq.[2]). 
The higher the SFR the more-massive is the upper limit of 
the ECMF. As massive clusters develop a top-heavy IMF 
according to eq. [3] a larger fraction of massive stars are 
formed with increasing SFRs. Fig. shows the dependence 
of cyigimf on the lower limit of the ECMF (M oc i,min)- For 
very flat ECMFs (dashed lines) this dependency is nearly 



6 C. Weidner, P. Kroupa, J. Pflamm-Altenburg 



Table 1. The IGIMF power-law slopes for a series of models with a changing 013 for clusters above Mi imit = 2xl0 5 Mq and for different 
lower limits of the ECMF, M cc i m ; n . The upper mass limit of the ECMF is chosen according to eq. [2] The slopes are obtained by a 
least-squares power-law fit to the IGIMF above 1.3 Mq. The columns designated as "dev" show the deviation of the top-heavy IGIMF 
from a pure power-law IMF. This is done by dividing the number of stars above 8 Mq from the top-heavy IGIMF model by the number 
for a pure power-law IMF with the fitted slope. Additionally, the columns entitled "ratio" list the ratio between the number of stars 
above 8 Mq divided by the total mass in stars. In the penultimate row, the mean SFR over a Hubble time (in Mq yr _1 ) is given, which 
would be needed to produce the total stellar mass over 14 Gyr instead of a burst of 10 Myr. In the final row the maximal cluster mass 
according to eq. [2] is shown, as well as the ratios between the number of stars above 8 Mq and the total mass in stars for = 2 and 
A^ecl.min = 5 Mq but without a top-heavy IMF in massive star clusters. 
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negligible. This is because for slops of ft less than 2.0 the 
ECMF is dominated by massive clusters and reducing the 
number of low-mass cluster by increasing A'/ cc i,mm does not 
affect the IGIMF significantly. For /3 = 2.35 the opposite is 
the case. For steep slopes, low-mass clusters dominate the 
ECMF and therefore have a strong influence on the IGIMF. 
Here, increasing M ec i jm j n changes the IGIMF strongly by 
eliminating the low-mass clusters. 



5 DISCUSSION 

With this contribution the SFR dependent IGIMF theory 
is extended into the starburst regime. The SFR dependent 
IGIMF, which produces results agreeing with observations 
for massive and dwarf galaxies |K roupa & Weidner 120031; 
Weidner fc Kroupal [20051 : iKoppen et al.1 120071 ; iRecchi et all 
2009t ). leads to top-heavy galaxy-wide IMFs in mas- 
sively starbursting galaxies when combined with the 
cluster-mass dependent high-mass IMF slope for glob- 
ular clusters and UCDs (eq. [3} by Marks, Kroupa 
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Figure 3. Dependence of the IGIMF slope above 1.3 Mq 
(chgimp) on the lower mass limit of the ECMF (M oc i im i n ). Solid 
lines show the results for (3 = 2.35 for SFRs from 10 (upper-most 
line) over 100 and 1000 to 10000 Mq yr" 1 (lowest line). The dot- 
ted lines are like the solid ones, but for /3 = 2.0 and the dashed 
lines for /3 = 1.6. 



& Dabringhausen (in preparation) and a SFR depen- 
dent lower limit of the ECMF or flat slopes of the 
ECM F. This is interesting be c ause some observational re- 



suits (|Finoeuenov et al. 


20031; Alonso-Herrero et al.l 2004 


Tornatore et alJ 120041: 


Loewensteinl 12006: Stanwav et al.l 


2005; Nayakshin & Sunyacv 2005: Alexander et al. 2007) in- 



deed suggest the IGIMF to be top-heavy when SFR > 100 
Mq j/r _1 . Current observational results from cosmological 
studies are very difficult to compare with the models because 
the observations do not measure IMF slopes and SFRs for 
individual galaxies but study indirect evidence (luminosi- 
ties, chemical enrichment, etc.) for whole populations and 
average the results over the galaxy luminosity function. The 
most-massive star cluster for which the IMF has been deter- 
mined by star counts, R136 in the LMC, sho ws no evidence 
for a non-Salpeter slope of the high-mass stars (|Selman et al.l 
1 19991 ). However, this cluster falls below the mass limit above 
which the parametrisation chosen here (eq. [3]) implies a top- 
heavy IMF. Additionally, no common description is used 
in the literature to characterize the top-heaviness. Some au- 
thors vary the peak of the IMF while others vary the slope or 
the lower and upper mass limit. Some available data points 
are shown in Fig. [4] together with the mode l predictions. 
The co smologically i nterpreted observat ions bv lvan Dokkuml 
l|2008l . triangle) and lBaugh et al. (2005, open circle) suggest 
rather flat IMFs already for relatively low SFRs and are dif- 
ficult to reconcile with the here presented model even for 
values of M ec i jm i n which are very high (10 s M e ). But as can 
also be seen in Fig. [4] the IGIMF results are in good agree- 
ment with the constraints for the Galactic and M31 bulge 
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Figure 4. IGIMF slopes above 1.3 Mq (qigimf) for the SFR val- 
ues of Tab. [T] and some observational constrains in dependence of 
the total SFR. The shaded region between the solid lines marks 
the range of model results for an ECMF slope f) = 2.0 while 
the dotted lines show the full en velope constr ained by all mod- 
els. The three asterisks are from iDavel ll2008ll . the one triangle 
around Q3 ss 0.7 is fr om Iva n Dokk uml ||2008|) . the filled circle 
with error bars is from lBallero et al.l J20 07all while the op en cir- 
cle with error bars corresponds to the Baugh et al. (2005) result. 
The dashed horizontal line marks the IWilkins et alJ {2008|) con- 
strain with the light dashed lines 0.15 dex above and below being 
their uncertainty range. 



{filled circle with er r or-bar , iBallero et al ] |2007al ). as well 
as the IWilkins" et alJ (^OOct ) constrains for the present-day 
mass density from the cosmological SFH (dashed lines). The 
GAMA-team finds a very similar trend of decreasing slope 
with SFR as the models in their sample of ~ 40000 galax- 
ies l|Gunawardhana et al.l[2010l ). Likewise, the IDavel (2008, 
asterisks) study of the amplitude of the galaxy stellar mass- 
star formation rate relationship seems to be in reasonable 
agreement with the model. Though for this agreement it is 
necessary that for galaxies with SFRs above 10 Mq yr" 1 
the ECMF has either a lower limit larger then 10000 Mq or 
a slope flatter then 1.5. Generally, it should be possible to 
test these scenarios observationally by observing the ECMF. 

Summarising, the main physical reason why the IGIMF 
becomes top-heavy at galaxy- wide SFRs > 10 Mq yr _1 is 
the formation of very massive star clusters with masses M cc i 
> 10 6 Mq. To achieve such high SFRs a galaxy with neu- 
tral gas needs to either globally be unstable (e.g. at early 
cosmological times when the gas fraction was very high) or 
be compressed globally due to an external tidal force field. 
Under these conditions the pressurised interstellar medium 
will collapse to massively giant molecular cloud complexes 
with masses of 10 s Mq or larger. In these complexes, ultra- 
compact-dwarf galaxy type star clusters may form with in- 
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dividual masses > 10 6 Mq, and these would be h aving top- 
heavy stellar IMFs ^Dabringhausen et al.l l 2009), as intro- 
duced into the IGIMF theory here. It may also be possible 
that at the same time the formed star-cluster mass func- 
tion becomes bottom light under extreme SFRs, which may 
be due to the conditions for forming low-mass (< 100- 1000 
Mq) clusters not being available, perhaps due to the intense 
stellar feedback which may suppress the formation of such 
low-mass molecular cloud cores. 

It should be noted he re that the two rather extreme 
top-h eavy IMF suggestions (|Baugh et al ]|2005l ; lvan Dokkuml 
2008) with ^empirical ^ 1 are difficult to reconcile with the 
existence of globular clusters today. To obtain such flat 
IGIMFs the massive star clusters would need to have such 
flat IMFs that more than 90% of their mass would be in 
stars which explode as supernovae. Such extreme mass-loss 
would quickly fully disperse these clusters, leaving no glob- 
ular clusters behind. 
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APPENDIX A: THE CANONICAL IMF 

The following multi-component power-law IMF is used 
throughout the paper: 

k ' (t^h") ™° ,m low ^m<m H , 

(^) ai ,m H ^ m < m , 

, fen fen ^<"<^ 

with exponents 



£(m) 



a = +0.30 
ai = +1.30 
a 2 = +2.35 
a 3 = +2.35 



0.01 < m/M Q < 0.08, 
0.08 sC m/Mg < 0.50, 
0.50 < m/M < 1.00, 
1.00 < iu/Mq < m max 



(A2) 



where dN = £(m)dm is the number of stars in the mass 
interval ratoml dm. The exponents Qj rep resent the stan- 
dard or canonical IMF (|Kroupall20ol . 1200a ). Though, a 3 is 
kept constant at 2.35 only for star clusters with M ec i less 
than 2xl0 5 Mq. For more massive clusters a 3 is changed 
with cluster mass according to eq. [3] 

The advantage of such a multi-part power-law descrip- 
tion are the easy integrability and, more importantly, that 
different parts of the IMF can be changed readily without af- 
fecting other parts. Note that this form is a two-part power- 
law in the stellar regime, and that brown dwarfs contribute 
about 4 per cent by mass only and that a discontinuity 
near mn impl ies brown dwarfs to be a se parate popula- 
tion (fc' ~ L iThies fe Kroupal [20071 . I2008D . A log-normal 
form in the stellar regime but below 1 Mq with a power- 
law ex tension to high masses was suggested by IChabrien 
(2003 ) but is virtually identica l to the canonical IMF (fig. 8 
in Dabringhausen et al.ll2008l ). 



The observed IMF is today understood to be 
an invariant Salpete r/Massey power-law slope (|Salpeterl 
1 19551 ; iMassevI 120031 ) above 0.5 M , being independent 
of the cluster density and metallicity for metallic- 

Isirianni et al.l 
120021, 



ities Z> 0.002 llMassev fc Hunted Il99 



2005 ; iWvse et al 



| 2000| . |2002|; iParker et al l l200ll; iMasseyj Il998l. 



2002; Bell ct al 



2003; iPiskunov et all 



l2004l ; |Pflamm-Altenburg fc Kroupall2006t ). Furthermore, un 
resolved multiple stars in the young star clusters are not 
able to mask a s ignifican tly different slope f or massive stars 
jMafz Apellanij|2008l ; I Weidner et al. l l2009l ). lKroupal (|2002l ) 
has shown that there are no trends with physical condi- 
tions and that measured high- mass slopes, 03, are a Gaus- 
sian distribution about the Salpeter value thus allowing 
us to assume for now that the stellar IMF is invariant 
and universal in each cluster. There is evidence of a maxi- 
mal mass for stars (m max * ~ 150 Mm, I Weidner fc Kroupal 
|2004|). a result later con f irmed by several indepe ndent stud- 
ies llOev fc Clarke! 120051; lFiger j l2005l ; lKoenll2006l ). However, 
according to Crowther et al. ( 2010l ) m ma x* may also be as 
high as 300 Mq. 
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